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Abstract 

We calculate the distribution of current fluctuations in two simple exclusion models. Although 
these models are classical, we recover even for small systems such as a simple or a double barrier, 
the same distibution of current as given by traditionnal formalisms for quantum mesoscopic con- 
ductors. Due to their simplicity, the full counting statistics in exclusion models can be reduced to 
the calculation of the largest eigenvalue of a matrix, the size of which is the number of internal 
configurations of the system. As examples, we derive the shot noise power and higher order statis- 
tics of current fluctuations (skewness, full counting statistics, ....) of various conductors, including 
multiple barriers, diffusive islands between tunnel barriers and diffusive media. A special attention 
is dedicated to the third cumulant, which experimental measurability has been demonstrated lately. 



1 Introduction 



A constant voltage difference across a conductor drives an electrical current which will always 
fluctuate around its mean value. Fluctuations result from random microscopic processes (thermal 
relaxation, scattering, tunneling...) undergone by the charge carriers. These fluctuations can be 
considered as an undesirable noise but also as a rich signature of the basic transport mechanisms 
occurring in the conductor. This second perspective has concentrated much attention in the 
mesoscopic community over the last decade [7]. 

In our previous paper [15] we gave evidence that the statistics of current fluctuations in a large 
classical model, the symmetric exclusion process, are identical to the ones derived for quantum 
mesoscopic conductors [27]. Here, we show that exclusion models allow also to recover the current 
fluctuations of small systems such as a single or a double barrier, even in the ballistic limit. 

In the present paper, we develop a classical approach to derive the statistics of current fluctua- 
tions in mesoscopic conductors ( "quantum conductors" ) and more generally in conductors smaller 
than the electronic inelastic mean free path and for some inelastic conductors. Solving the current 
statistics problem is reduced to finding the largest eigenvalue of a modified evolution matrix, later 
called the counting matrix. We recover and generalize the well known current statistics for a few 
mesoscopic systems. Our description is based on the exclusion process models, which have been 
widely studied in statistical physics and probability theory [46, 17, 30]. The main benefits of this 
approach are its conceptual and analytical simplicity. 



In the remaining part of this introduction section, we briefly recall the traditional approaches 
for mesoscopic transport (section 1.1) and the basic mathematical tools necessary to describe 
current fluctuations (1.2). Section 2 presents two exclusion models fitted for condensed matter 
conductors and the procedure to derive the complete statistics of current fluctuations (later called 
"Full Counting Statistics" or FCS.). In section 3, our exclusion models are used to derive the 
current statistics of various elementary conductors. 



1.1 Traditional formalisms for transport in condensed matter physics 

A number of approaches have already been used to describe the FCS in mesoscopic conductors. 
The Scattering Matrix theory[8, 43, 9, 28, 32, 27] is well adapted to the modeling of quantum- 
mcchanically coherent conductors in a regime where electron interaction effects are sufficiently 
weak to be neglected. With this strong assumption, this allows to treat an arbitrary large number 
of transverse conduction channels, which yield independent contributions to the current statistics. 
This approch, combined with results from random matrix theory for scattering matrices [3] has lead 
to precise predictions for the FCS of a disordered conductor in the diffusive regimc[25]. A more 
direct microscopic treatment of disordered systems relies on the Keldysh technique [21] to construct 
the non-equilibrium density matrix of the steady state at finite current. Disorder averaging is then 
performed using a non-linear Sigma model representation[20]. A rather general circuit theory has 
been constructed to account for the influence of an arbitrary environment, described in terms 
of an equivalent circuit, on the measured fluctuations of a mesoscopic conductor [36, 38]. Semi- 
classical descriptions, based on the Boltzmann-Langevin model[44, 18] have also been used to 
derive the first four cumulants of current fluctuations in a diffusive medium [33, 34]. Other semi- 
classical approaches focused on high order statistics and FCS of a double tunnel barrier [12], chaotic 
cavities[39] and diffusive media[42, 15]. The semi-classical results are the same as the ones obtained 
with the corresponding quantum conductor model. 

The exclusion models discussed in this paper represent an extreme semi-classical approxima- 
tion : the only quantum rule which is preserved is Pauli exclusion principle. In particular, electrons 
have no phase and don't interfere. 



1.2 Mathematical formalism for Current fluctuations 

If q t is the algebraic charge which flows accross a section during time t, the fluctuations of current 
I = q t /t depends in principle on the duration t chosen to measure /. In practice the long time 
response of the measuring electronics apparatus sets a lower bound on t : this bound is most often 
decades larger than all the physical times experienced by charge carriers (diffusion time, dwell 
time, coherence times in the conductor and in the electrodes,...). Thus, experiments correspond 
to the t — > oo limit, often called the zero-frequency limit in the shot noise literature. In this limit, 
the choice of the cross-section is irrelevant since the maximum charge accumulation between two 
different cross-sections is finite, at least in a conductor connected to two electrodes only. 

In a conductor smaller than the inelastic mean free path, carriers do not undergo inelastic 
collisions. It is therefore reasonable in many situations to neglect interaction effects on such small 
length scales. Equivalently, we may then assume that these charge carriers remain on independent 
energy levels [7]. Consequently the statistics of the total current will consist in a summation of 
independent random variables corresponding to different energy levels. In the following, to keep 
equations free of elementary-charge prefactors, we focus on carriers counting rather than charge 
counting. In addition, we will call this generic charge carrier an electron. 

If Pt,t{Q) is the probability that Q electrons have been transfcrcd at the energy level e during 
a time interval t, one can fully characterize the counting statistics by the cumulant generating 
function : 



StA z ) = ln 

or equivalently by cumulants (the n th order one is written here) : 

d n S t A< 
d(lnz) n 



In(zQ) (1) 



= (2) 



we have in particular : d = Q , C 2 = (Q - Q) 2 , C 3 = (Q - Q) 3 , C 4 (e) = [Q - Q) 4 - 

For a given conductor, the current at an energy level e only depends on the boundary conditions, 
that is the fillings p L {e) and p R (e) of the left and right electrodes (or "reservoirs") at both ends of 



the conductor. If we rewrite the cumulant explicitly as C n (t,e, p L , p R ), the cumulants K n (t) for 
the whole conductor are given by 



K n (t) = J C n (t,e,p L (e),p R (e)) n(e) de (3) 

where n(e) is the density of energy levels in the conductor. Likewise, the cumulant generating 
function for the whole conductor can be derived with the same type of summation. For comparison 
with experiments Fermi-Dirac distributions are imposed in the left and right electrodes : 

P L (e) = l ~7^v ; P R (e) = : (4) 

l + e i^T l + e k BT 

With such fillings, the K n (t) are function of the driving voltage normalized by temperature 
eV/ksT. The fc^T » eV limit corresponds to the Johnson-Nyquist thermal noise and the op- 
posite limit to pure shot noise. In this paper, the integration Eq. (3) over e will be estimated 
assuming that n(e) and C n {t 1 e, p L , p R ) are independent of e. This assumption is quite reasonable, 
since in most cases the Fermi energy in the reservoirs is much larger than both the thermal energy 
window fceT and the driving energy eV. 

The electrical conductance G and the current noise power density Si are proportional to K\ 
and K 2 : 

G = I/V = eKi/Vt ; S 1 / = 2 J SI(T)SI(0)d,T = 2e 2 K 2 /t (5) 

where SI(t) = I(t) — I is the current fluctuation at time r. The time scale of the model dynamics 
can be chosen arbitrarily since this only changes the prefactor of the cumulant generating function. 
The transport mechanism is characterized by the cumulants C2, C3,... (or K 2 , K^,...) normalized 
by C\ (or K\). In particular, we will focus on the normalized shot noise power 

F = 5j/2eT = Ki/Ki (6) 

and the normalized skewness 

F 3 = Ka/Id (7) 



2 Exclusion Models 



In this section, we first present an exclusion model mostly adapted to the modeling of nearly 
ballistic conductors. We call it the counter-flows exclusion model because the two directions of 
propagation of electrons found in a ID conductor are explicitly considered. Conductors with a 
low transmission efficiency, such as tunnel barriers or diffusive media, can sometimes be described 
by a simpler exclusion model, presented in the tunnel exclusion model section. Many systems 
studied in the exclusion, hopping-model and sequential-tunneling literatures are directly relevant 
to this latter category of conductors. These models describe independent particles, apart for the 
exclusion constraint which represents the effect of the Pauli principle. 



2.1 The Counter-Flows Exclusion model 

The counter-flows model is inspired from the Landauer[26, 48, 8, 24] picture of conductors : at 
zero temperature, electrons are injected periodically from the reservoirs to the conductor. This 
assumption seems a good enough modeling to account for the FCS in the t — > 00 limit. Indeed, 

the predictions of the model would remain unchanged if the variance (N t — N t ) 2 of the number 
of injection attempts N t during a time interval t, is only sublincar in t. This later property fol- 
lows from the Pauli exclusion in degenerate electrodes which imposes an anti-correlation between 
injection events[28, 27]. The experimental validation of Landauer approach[41, 23, 7] justifies a 
posteriori this nearly-periodic injection model. While in the sample, these charge carriers may 
undergo internal scattering on localized barriers and finally are either reflected or transmitted to 
electrodes at both ends of the conductors. The Pauli exclusion principle is fulfilled at each stage 
during the system evolution. In the original model, electrons are described by wave-packets and 
phase coherence is preserved during the scattering but alternative versions of this model dropped 
the phase information (for example, see[13, 31]). 

More precisely, the 27V sites counter-flows model consists in N + 1 barriers, each characterized 
by 2 transmission probabilities (from left to right) and (from right to left) where i is 
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Figure 1: Counter- flows model. Probability of evolution for the various configurations of electrons 
reaching the i th barrier at time t. 
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Figure 2: Counter-flows model. Upper Fig. : Transmission and reflection probabilities for the i 
barrier, assuming no conflict with the exclusion principle. Lower Fig. : The counter-flows model 
for N = 4. The white circles represent empty sites, the black disks are electrons, the gray disks 
stand for sites with a fixed filling probability and the arrows indicate the direction of propagation 
associated with each site. 



the index of the barrier increasing from left to right (1 < i < N + 1). Between two consecutive 
barriers, 2 sites are available for at most 2 electrons propagating in opposite directions. So a 

configuration at time t is characterized by 27V binary variables (t) and ■* (t) for 1 < i < N; 

Tj (t) (respectively T^'(t)) is equal to 1 if an electron propagating to the right (respectively 
to the left) is present at site i at time t. Time is discrete and at each time step, electrons are 
transmitted through one barrier to the next site, unless a back-scattering occurs on the barrier. 

By definition of the dynamics of the model, rfl^(f + 1) and T^\t + 1) depend only on T^\t) 

and on -rfc (t), and the (classical) transition probabilities are given in figure 1. This allows for a 
simultaneous update of all occupancies, even in the presence of backscattering on barriers. 

At the boundaries of the conductors, each electrode is modeled by 2 sites, the occupation 
states of which are re-set before each time step. The site corresponding to an electron propagating 
into the conductor is re-filled with probability p L (left electrode) or p R (right electrode) and the 
site accessible to the electron leaving the conductor is re-emptied at each time step (see Fig. 2). 
The densities p L and p R are given by Fermi-Dirac distributions (Eq. (4)). After this reset, the 
one-time-step evolution follows the same transmission/back-scattering rule that holds in the bulk 
of the conductor. 

On modeling real conductors, the barriers can represent junctions (between two different mate- 
rials for example), scattering centers (impurities, structural defects, ...) or even inelastic processes 
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Figure 3: Tunnel model. Upper Fig. : Tunneling probabilities across the i and (i + l) th barriers 
for an electron located between them, assuming no conflict with the exclusion principle. Lower 
Fig. : Tunnel exclusion model for N=5 sites. 
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Figure 4: SSEP model for N=5 sites. T L , T and T R are the tunneling probabilities and p L , p R the 
electrodes' fixed fillings . 

(phonon or photon- assisted hopping, emission of phonon or photon, ...). The model parameters N, 

T^~^ and ^ are related to the corresponding physical quantities such as tunneling probabilities 
or scattering cross-sections. 



2.2 The Tunnel Exclusion Model 

It is useful to note that the counter-flows exclusion model may be decomposed into two indepen- 
dent stochastic models. Let us define new variables <Ji(t) and o-[(t) such that <7j(t) = T^\t) if 
% and t have the same parity and else o~i(t) — (t). In a similar way, cr^(t) = t\ (t) if i and 
t have the same parity and else <j[(t) = t\ (t). From the definition on the model, the random 
variables <jj are completely decoupled from the a[ variables. It turns out that in the limit of small 
transmission probabilities, the dynamics of each of these two ensembles of binary variables may 
be formulated in terms of a simpler lattice model, that we shall call the tunnel exclusion model. 
The elementary time-step of the latter model involves two steps in the former one. This has the 
advantage that in the limit of a vanishing transmission probability, the configuration of oVs does 
not evolve in time. For each of these two independent submodels, expanding the evolution of this 
reduced system to first order in transmission probabilities, and taking the continuous time limit, 

we get the model which definition is sketched on Fig. 3. In this case, the quantities and 
become the probabilities per time unit of tunneling across the i th barrier from left to right and 
vice-versa, provided that the target site is empty. Each electrode is modeled by a single site, the 
occupation of which is reset to p L (left electrode) or p a (right electrode) before each time step. 
The fillings p L and p R are given by Fermi-Dirac distributions (See Eq. (4)). 



A special choice of the tunneling probabilities is the Symmetric Simple Exclusion Process 
or SSEP [45] (see Fig. 4) for which the the internal barriers are symmetric (t[^ — and 
uniform along the conductor (independent of i for 2 < i < N). We note this probability T. 
The two out-most barriers are also modeled with symmetric rates T L — T± and r s = rjv+i. 
Physically they account for the electrical connection between the electrodes and the conductor. 
In the theory of exclusion processes [15], one usually represents the reservoirs by injection rates a, 
S, and extraction rates 7, (3 which give an equivalent description of the boundary conditions if: 



a = p L r L ,5 = p R r R , 1 =(l-p L )T L ,(3=(l-p R )T R . 



2.3 The FCS solving procedure 



N 



In the counter-flows model, the conductor has 2 2N internal configurations C = {r 1 ( ^ ) , ... 

Let Pt(C) be the probability of finding the system in configuration C at time t. As the dynamics 
is a Markov process, the evolution equation for pt(C) can be written : 



p t+1 (C) =J2[Mi(C,C') + M (C,C')+M_ 1 (C,C')}p t (C) 



(8) 



where we have decomposed the evolution matrix into three parts Mi, Mq and M_i, depending 
on whether, when the system jumps from configuration C to configuration C, the total number of 
transfered charges increases by 1,0 or —1. 

If we define Pt(C, Q) as the probability that the system is in configuration C at time t and that 
Q charges have been transfered, one has : 

P t+ i(C, Q) = ]T Mi(C,C')P t (C, Q-l) + M Q (C,C')P t (C, Q) + M_i(C,C')P t (C , Q + l) (9) 
c 

Then the generating functions Vt(C, z) defined by : 



satisfies 



V t+ i(C,z) = J2 



V t (C,z)= ]T P t (C,Q)zQ 

Q=-oo 



z Mi(C,C) + M (C,C) + -M_i(C,C) 



V t (C',z) 



If we introduce M z that we will call the counting matrix, defined by: 



(10) 



(11) 



M Z (C,C) = z Mi(C,C) + M {C,C) + -M_i(C,C) 



(12) 



it is clear from Eq. (11) that in the long time limit, the cumulant generating function for the total 
number of transfered charges is : 



S t (z)=ln(zQ)=ln 



ln(u(zY) ~t\n(v{z)) 



(13) 



where v{z) is the largest eigenvalue of the counting matrix M z [16, 15]. Due to the fact (see 
beginning of section 2.2) that the counter-flows model can be decomposed into two decoupled sets 
of variables, the eigenvalue v(z) can in fact be obtained by diagonalizing a 2 N x2 N matrix. 

In the tunnel model, the conductor has 2^ internal configurations and -as previously- we 
call pt{C) the probability of finding the system in configuration C at time t. The time being 
continuous in this model, one has : 



dp t (C) 
dt 



- Y}Wi(C,C) + W (C,C) + W^(C,C')} Pt (C') 



(14) 



where the evolution matrix has been decomposed into three parts Wi, Wq and W-i, depending 
on whether when the system jumps from configuration C to configuration C, the total number of 
transfered charges increases by 1,0 or —1. Eq. (14) is a continuous time version of Eq. (9), the 
main difference being the diagonal elements of Wq are now all negative. 

Following the same procedure as above, we can define the counting matrix W z by : 



W Z (C,C) = z Wi(C,C) + W (C,C) + -W-i(C,C) 



(15) 



and we find the cumulant generating function for the total transfered charge in the long time 
limit : 

(16) 



S t (z) = ln(zfl) - In (e^ *) - t p,{z) 
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Table 1: Asymmetric Single Barrier for arbitrary fillings p L and p R of the electrodes. 



where p(z) is the largest eigenvalue of the counting matrix W z . This latter equation can be seen as 
the first term in the expansion of the corresponding equation obtained in a discrete time approach. 

Both for the counter-flows and tunnel models, the FCS is fully determined by the largest 
eigenvalue of what we called the counting matrix. 

The full knowledge of the eigenvalue is not necessary if only the first n cumulants are wanted. 
In this case, the equation satisfied by the eigenvalues \M Z — v{z)T\ = (or \ W Z — p,(z)l\ = 0) can 
be solved by a perturbation theory and the n th cumulant is obtained from the coefficient of the 
Q th or der of the eigenvalue in powers of log(z) (see Eq. (2)). Once the counting matrix is written 
down, this procedure can be easily performed by an analytical calculation software. 



3 Application to mesoscopic systems 



In the remaining of this paper we derive the FCS or the first cumulants of basic mesoscopic sys- 
tems. A special attention is dedicated to the current fluctuations skewness (third cumulant), the 
physical intcrcst([29, 19]) and measur ability [40] of which have been recently emphasized. Indeed, 
at high temperature the skewness can reveal information about transport which are not blurred 
by thermal fluctuations[29, 19]. Some of the results derived are already known and they validate 
exclusion modeling for charge conduction in condensed-matter systems. The various new results, 
often derived in a few lines of linear algebra, illustrate the strength of this modeling. 



3.1 Asymmetric barrier and Single channel 

The counter-flows model with N = site is a single barrier between two electrodes of fillings p L 
and p R . Since the system has no internal state, the counting matrix M z is a scalar. A positive 

charge transfer (from left to right) will occur with probability p + = p L (1 — p R )T^^ and a negative 

transfer with probability p_ = (1 — pi,)p R T^*~^ (see Table 1). 

Following the general procedure for the counter-flows model, we consider the counting "ma- 
trix" : 

M z = P+ z + p_ z- 1 + (l- P+ -p_) (17) 



The logarithm of the largest (and unique) eigenvalue v(z) = M z gives the cumulant generating 
function S t (z), which fully characterize the FCS of an asymmetrical barrier : 

s t {z)/t = \n{[p L (i-pjrH] (z-i)+ [(i-pjp«r<^>] O^-rj+i} (is) 



A few particular cases are interesting : 

• The ]?( <_ ' — > 1 limits account for quasi-ballistic barriers. On the opposite case of tunnel 
barriers (r^T^ < 1), the FCS can be re-estimated from a first order expansion of Eq. (18) 
or directly with the tunnel exclusion model with N = : 

S t (z)/t = n{z) = W z = p L (l- p R )T^\z - 1) + (1 - pjp/^'fz- 1 - 1) (19) 



• The asymmetric case (T^~ ^ =/= T^~^) accounts for inelastic barrier, such as those for which 
stepping over the barrier requires the emission or assistance of a photon or phonon[29]. 



For symmetric barriers 

r = = r<«-> (20) 

we recover the important case of a conduction channel of transparency T encountered in meso- 
scopic transport [28]. This is to be expected from our discussion in section 4, where we show 
that for a single barrier, the counter-flows model can be precisely related to the evolution of the 
diagonal part of the density matrix of a quantum many-electron system. Once the behavior of 
a single conduction channel has been determined, scattering matrix theory shows how to reduce 
the problem of interaction-less electronic transport through a quantum constriction into a set of 
independent symmetric barriers. In the zero temperature limit, the only states which contribute 
to the FCS are those whose energy e is such that p L (e) = 1 and p R (e) = 0, or p L (e) = and 
p R (e) = 1. The FCS is then given as a superposition of independent binomial laws ("partition 
noise"), one for each of these scattering states, lying in an energy window eV, where V is the 
voltage drop accross the barrier [28]. In the high temperature limit, one has to integrate the single 
channel result Eq. (18) over the complete Fermi-Dirac distributions in the reservoirs, given in 
Eq. (4). 

From equations (2), (3), (18) and (T = r^ - ^ = r^)), one can derive the normalized noise 
power F — K2/K1 in the low temperature limit and the normalized skewness F3 = K3/K1 in 
the high temperature limit [29, 20]. It is interesting to note that these two quantities turn out to 
be equal. More generally, for a mesoscopic conductor decomposed into independent channels of 
transparencies % one has : 

F(eV » k B T) = F 3 (eV « k B T) = ^^T^ (21) 

The physical information contained in the third cumulant at high temperature[29, 19] is the same 
at the one contained in the low temperature second cumulant. 

Eq. (21) will be directly checked for the mesoscopic systems considered in the rest of this paper. 



3.2 Double barriers 



For single barriers, the agreement between the exclusion model and mesoscopic models is not 
surprising since the boundary conditions (injection from the electrodes,...) are identical. For 
double barriers, it is crucial to account properly both for the boundary conditions and for Pauli 
exclusion principle inside the conductor. 

Double barriers have been widely studied because a rich behavior results from the interplay of 
various effects including Pauli exclusion principle, Coulomb interactions, inelastic processes and 
quantum resonance [7]. In this section, we first present some results on a generic double barrier. 
Then we see how this system relates to various experimental devices (quantum dots, hopping on 
localized states, islands and wells) and how the Coulomb interaction between electrons can be 
introduced to account for charging effects. The case of a partly or fully diffusive island between 
two tunnel barriers is addressed in section 3.3. 

Generic double barrier 

We first consider two symmetric barriers of transmission Ti and T 2 , temporarily in the zero 
temperature limit eV ^> k B T (p L = 1 and p R = 0). The upper graph of Fig. 5 depicts the 
corresponding counter- flows model, with N = 1, while the lower graph labels the 2 2N internal 
states of the system. 

If the charge counting is done over the second barrier (arbitrary choice), the counting matrix 
M z is : 

/ T^z T 1 Y 2 z 1 \ 

(i-ri)r 2 * i-Tx 
ri(i-r 2 ) i-r 2 
\ (i-r 1 )(i-r 2 ) j 



M z = 



(22) 



In this matrix, the states are ordered from state- A (upper- left) to state-D (lower-right). The 
eigenvalues of M z can be easily found and the cumulant generating function St(z) is proportional 
to the logarithm of the largest one : 



S t {z)/t = \n 1- 



ri + r 2 -r 1 r 2 ^ 




1 - 



ri + r 2 -r 1 r 2 z 



(i-ri)(i 



(23) 
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Figure 5: Upper Fig. : Double symmetrical barriers (counter-flows model with AT = 1, = rj - ^ = 
r,). For legibility, the reflection probabilities 1 — T\ and 1 — T 2 are not written. Lower Fig. : The 
internal states of the system. 

The symmetry of this expression between Ti and r 2 illustrates that the charge counting can be 
performed on any side of the system without changing the result. As expected, the single barrier 
FCS is recovered if one barrier is transparent (FiorF2 = 1). Eq. (23) extends two other results 
first obtained by de Jong in the Boltzmann-Langcvin formalism : the first and second cumulants 
of a double barrier[13] and the FCS of a double tunnel barrier (rir 2 -C 1) [12]. It would be 
interesting to relate these classical expressions to those of a full quantum treatment of the double 
barrier system. As described below, such a precise connection may be established in the tunneling 
regime where both transmissions T% and r 2 are very small. 

For arbitrary fillings p L and p R of the electrodes, the counting matrix M z has no zero element 
and the eigenvalue problem is still manageable but more tedious. The power expansion method 
presented in section 2.3 is chosen to derivate the cumulants C\, C 2 and C3, and integration 
over Fermi-Dirac distribution in the electrodes, according to Eq 3 and Eq 4, gives the cumulants 
K^eV/kuT), K 2 {eV/k B T) and K 3 (eV/k B T). It is useful to define: 



One finds 



ri + ra-r^ (24) 



Cx/t = {p L -p R )^ (25) 

A 12 



n u + Pr ) - (pi + P l)TjTt(2 - r 12 ) - 2p L p R r 1 T 2 (r 2 1 + rj - r } r 2 (Xi + r 2 )) 

1 12 

(26) 

For eV <C k B T, we find K2/K1 — 2k B T/eV, that is the well known Johnson- Nyquist thermal 
noise formula, more often written Sj — Ak B TG where the conductance G and the current noise 
power spectral density Si are given by Eq. (5). 

We focus now on the normalized skewness F 3 = K3/K1. We give below its eV ^> k B T 
and eV <C k B T limits, which -as would be expected from a quantum mechanical derivation 
(Eq. (21))- is equal to the normalized noise power F = K2/K1 in the zero temperature limit. 
Fig. 6 shows that depending of ri and r 2 , the skewness changes sign in the low temperature 
limit (left Fig.) and not in the high temperature one (right Fig.). The insert on the left Fig. 
shows F 3 = (1 — r 1 r 2 /r 12 )(l — 2r 1 r 2 /r 12 ) obtained when the exclusion principle is deactivated 
between the two barriers. The change of sign is still observed and thus, it cannot be attributed to 
correlation inducted by the exclusion principle inside the conductor. 



F 3 (ev » k B T) = r f + r l + r ? r i( 4 + r i + r 2) + £gr|(6 - zv t -T\- 3r 2 - 1%) - r 1 r 2 (2r? + rf + 2r| + ri) 



r 4 
1 12 



(27) 

F 3 (eV « k B T) = 1 - ElI^LzIlH] = F{eV » fcsT) (28) 

1 12 



Tunnel limit 



The FCS of a double tunnel barrier can be obtained from the general expression Eq.(23) for 
the double barrier. Nevertheless, it is interesting to derive it with the Symmetric Simple Exclusion 




Figure 6: Normalized skewness F3 = K3/K1 for a double symmetrical barrier. Left Fig.: Zero 
temperature limit [eV 3> ksT) (insert: the exclusion principle is deactivated between the barriers). 
Right Fig.: High temperature limit (ksT 3> eV). 

Model with N = 1. For arbitrary fillings p L and p R of the electrodes, and charge counting done 
over the second barrier, the counting matrix is : 

W z =( r r f~ P , R) l T vn~ pL \ TlP r L+T2P r Z ") (29) 
\T 2 (1- p R ) z + Ti(l- p L ) -T 1 p L -T 2 p R J 

Following section 2.3, the cumulant generating function S t {z) is proportional to the largest eigen- 
value of W z ■ After a few lines of algebra, we find : 




St{z)/t = —4-^ + \l ( -4r^ ) + { Pl (i - Pr ) (z - 1) + Pr (i - P J (z- 1 - 1)) 

(30) 

We recover the FCS derived in [2], which slightly differs from the one derived in [5]. 
Experimental systems 

Depending on the elastic or inelastic nature of the barriers, several double barrier systems are 
traditionally considered. If we restrict ourselves to elastic barriers, the size of the island between 
the barriers is another source of experimental diversity. 

In large but finite-size islands ("quantum island"), the quasi energy levels are discrete and 
the generic system described above can model each such level. At this stage, it may be useful 
to describe in more detail the connection between a microscopic quantum coherent model for 
a single conducting channel in the presence of a double barrier in the tunneling limit, and the 
corresponding SSEP. At zero temperature, the generating function St(z)/t for the microscopic 
quantum model is given by [28] : 



S t (z)/t = v F ^-ln(T(fc)(z-l) + l), 



where T{k) is the transmission coefficient for an incoming electronic wave-function with wave- 
vector k and fc max — fc m m = vf being the Fermi velocity in the electrodes. For a fully 

coherent system, T{k) has to obey the quantum series composition rule for transmission matrices. 
In the limit of small transmissions T\ and 72, this yields the well-known Lorentzian resonance 
profile: 

T (k) = ATl% 
K ' (T 1 +r 2 ) 2 + 16(fc-A: ) 2 Z 2 ' 
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Figure 7: Left Fig.: Energy levels associated with the two sites in interaction. Right Fig.: The 
three internal states of the system. The black disks are electrons and the circles are available sites. 



where I is the distance between the two barriers. In the limit where eV/h is much larger than the 
decay rate 7 = ^(Ti + T-i) of the resonant level in the central region, and if the average energy 
of this level is such that k falls in the interval [k min , fc max ] , we may take the k integration over 
the whole real axis, yielding exactly Eq (30) in the special case p L (e) = 1 and p R (e) = 0, provided 
we set Ti — ^% (i — 1, 2) [12] . So in the tunnel limit, the simple SSEP may be viewed as the 
result of integrating the contributions of quantum-mechanically coherent electrons over an energy 
window larger than the width fvy of a discrete level inside the cavity delimited by the two barriers. 

When the island is infinite is the transverse direction ("Quantum well"), the energy levels 
become a continuous energy band : the generic model still applies but the integration Eq. (3) 
should be done versus wavevectors[7] k z (longitudinal) and kj_ (transverse direction) rather than 
versus the energy e. 

Finally, when the island is small ("Quantum dot") or when it consists in localized states 
(dopants, impurities...), the energy levels are well separated but one can no longer neglect that 
electrons entering or leaving the island are changing its Coulomb electrostatic energy, which induces 
a shift of the energy levels ("charging effect"). We illustrate in the simple example below how to 
account for such an effect. 

Charging effect 

We consider now a localized state tunnel-coupled to two electrodes with Fermi-Dirac fillings 
p L (e) and p R {t) (Fig. 7). The coupling is elastic (for example, it could be due to some overlapping 
of the localized state wavefunction with electrodes) and characterized by two symmetrical trans- 
mission coefficients ri = T2 = T with no dependence with e. Up to two electrons can be trapped 
on the localized state (spins up and down for example). If no electron is trapped, the two available 
sites are degenerated at energy level eo. The trapping of one electron shifts to e\ = e + e 2 /C the 
energy level of the remaining empty site, where C is the effective capacitance of the island. This 
system has three internal states, corresponding to 0, 1 or 2 trapped electrons (States A, B and C 
on Fig. 7) and the corresponding counting matrix is 

/ -2p L (e )-2p R (e ) (1 - p h (e )) + (1 - p R (e ))z 

W Z =T\ 2p L (e ) + 2p R (e )/z -2 + p L (e ) - p L (ex) + p R (e ) - p R {e x ) 2(1 - p h (e x )) + 2(1 - p R (e x ))z 
\ p L (e x )+p R (e x )/z -2 + p L {e x ) + p R (e x )) 

(31) 

With a software like Mathematica, the analytical equations for the cumulants Ci, C2 and C3 are 
easily derived by the series expansion method (section 2.3). These expresssion are not reported 
here but have been used to generate the normalized cumulants of Fig. 8 (Thick lines) for parameters 
eo/ksT = 10, (e 2 /C)/fcsT = 20 and T <C 1. The first cumulant C x is normalized by the averaged 
transfer charge 2.iT/2 expected in the high driving voltage limit. For comparison, alternative 
models are shown on the same figure. The dashed lines correspond to 2 sites at fixed energy 
levels eo an d ei and connected to the electrodes through tunnel barriers of transmission T. This 
modeling of the original problem is expected to account for the high eV limit (ei <C eV). The thin 
line corresponds to 1 site between tunnel barriers of transmissions 2r (left barrier) and T (right 
barrier). This model is expected to account for the intermediate plateau regime for which only one 
site of the localized state contributes efficiently to transport (eo <C eV <§; ei). Fig. 8 shows that 
the general solution to the first problem (thick line) is in good agreement with the two limiting 
cases considered for comparison. 



3.3 Multiple barriers : from double wells to diffusive islands between tun- 
nel contacts 

In this section we discuss the general SSEP model presented at the end of section 2.2 , (N — 1) 
identical barriers of transmission f< 1 are sandwiched between two tunnel barriers of transmission 




Figure 8: Two sites interacting via the Coulomb electrostatic repulsion and tunnel coupled to 
electrodes (see text). Left fig.: Normalized current Ci/tT. Central fig.: Normalized noise power 
CijC\. Right fig.: Normalized current skewness C3/C1. Thick line : present model, thin and dash 
lines : alternative models for comparison (see text). 



T L and T R . Writing down that the mean current and its variance do not depend on the barrier 
over which they are measured, we showed in[15] how to derive the first two cumulants : 

Ci/t = ^(p L -p,) (32) 



2 / 1 

(Pl + Pr - 2 PlPr) - tAPl - Pnf 1 - 



3 vri - rR/ V. 2N i 2iVi 2 (iVi-l) 



(33) 



with Nl =N-l + ±- + ±- and A=£(£-l) ( 2 £ - l) + £ (£ - l) (2^ - l) 

For arbitrary Ni, the third cumulant C3 can be derived from the third order series expansion 
of the cumulant generating function given in[15] (Eq. (C.14)). We give below the normalized 
skewness in the large N\ (or TV) limit. 

Experimental systems 

Three limiting cases of the above general formula could be relevant to mesophysics. 

• Firstly, for N=2, the system is a triple barrier with three different tunnel transmissions. Such 
systems have been widely studied as a model for double localized state, coupled quantum dot and 
double well structures (for example see [7, 22]). It would be interesting to compare our predictions 
to the experimental results. 

• In the limit N — > 00 for constant T, T L and r R , the relative contribution of the two out- 
most barriers is vanishingly small and we obtain a diffusive medium in good contact with the 
electrodes. Section 3.4 addresses this purely diffusive regime and gives its FCS. At this point, we 
just point that C\ and C2 are in agreement with the one found by alternative condensed-matter 
formalisms [33, 4, 35, 1, 13, 6, 47]. 

• The limit N — > 00 for constant T, iVT^ and NT R accounts for a diffusive island between 
two tunnel contact, the resistances of which are comparable to the resistance of the diffusive island 
itself. We define as q L = T/(NT L ) and q R = T/(NT R ) the ratios of the left and right contact 
resistances over the diffusive island resistance. The large q L and q R limit accounts for simple or 
double tunnel barriers while vanishing values of q L and q R account for a purely diffusive medium. 

The normalized noise power F (Eq. 6), is found to be, for an arbitrary temperature : 



with 

3 3 (1 + q L + q 



a = t+ _ ^" ^ n3 (35) 
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Normalized cumulants : (tr/ > N) 
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Table 2: First cumulants C n for diffusive medium for arbitrary filling of the electrodes, and low 
and high temperature normalized cumulants K n /K\. (X = (p L — p R ) and Y = (p L + p R — 1)). 



The well known noise power of double tunnel barriers[10, 11, 7] and diffusive media[33, 4, 35, 1, 
13, 6, 47] are recovered in the large and small q L , q R limits. The normalized skewness is found to 
be : 



F, 



K3 
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(36) 



(37) 



Eq. (36) bridges continuously between two simple situations : double tunnel barriers (large q L , 
q R ) and diffusive media (small q L , q R ) for which we also recover known results[20, 34]. In the low 
and high temperature limits, we found : 



F 3 (eV^>k B T) = l + 2f 



F 3 (eV « k B T) 



F{eV > k B T) 



(38) 



The high temperature skewness is equal to the low temperature noise power, as would be expected 
from a fully quantum treatment (see Eq. (21)). 



3.4 Diffusive medium 



In this section, we consider a diffusive medium with good contacts to the electrodes. We first 
consider the SSEP model in the N — > oo limit, which is enough to account for a purely diffusive 
behavior. At the end of this section, we consider the cross-over from a ballistic to diffusive 
conduction. To do so, the counter-flows model will be required. 

SSEP diffusive medium 

In our previous paper[15], we detail the derivation of the FCS for SSEP models in the N — > oo 
limit and we shall only summarize the derivation here. The cumulant generating function St 
depends on N, z, pl and pr but to first order in 1/N, S t turns out to depend only on a combinaison 
lo of these parameters : 

w = (f - ^){PLZ - Pr- PLPnjz - 1)) ^ 



Consequently, to first order in 1/N one can write 

S t ft 



(40) 



where the expression of R(u>) has been conjectured. Thus, if the FCS is known for two arbitrary 
fillings of the electrodes, then it can be deduced for any fillings or at any temperature. 

Based on the exact derivation of the first 4 cumulants, we conjectured that for p L = p R = 1/2, 
the statistics of current fluctuations is Gaussian (at order 1/N). This fully determines R(w) and 
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Figure 9: Normalized current skewness K3/K1 in a ballistic/diffusive medium parametrized by 
the ratio L/L e i ast of the conductor's length L over a elastic mean free path L e i ast (counter-flows 
with N = 8). The two lines correspond to known limit for L/L e i ast 3> 1 (continuous line) and for 
L/ Leiast <C 1 with eV 3> ksT (pearly line). 

thus the current statistics at a arbitrary fillings p L and p R : 

s t /t = ^[io g (VTT^ + V^)] 2 (41) 

This expression is identical to the one found at zero temperature by[25] and at arbitrary 
temperature by[20] (and by[37] with minor discrepancies). Tables 2 gives the first cumulants C n , 
K n /Ki for ksT <C eV and for ksT ^> eV. Right/Left symmetry and particle/hole symmetry 
suggest the change of parameters : X — (p L — p R ) and Y = (p L + p R — 1). Agreement is found 
with the first four cumulants which can be easily derived from[34]. 

Ballistic-Diffusive cross-over 

We now consider a counter-flows model with uniform symmetrical transmission coefficients T 

in the bulk of the conductor (T = r\~*^ = T^^ for 1 < i < N + 1 ), unity transmission at 

the boundaries = = rfc^ = rfo^ = 1), and we consider the N — > 00 limit taken 

at constant (1 — T)(iV — 1). We can define the mean free path L e i ast of a single electron as 
1/(1 - r) and the length L of the conductor as N - 1. Then (1 - T)(N - 1) = L/L elast < 1 
corresponds to a ballistic conductor while L/L e i ast ^> 1 to a diffusive one[14, 31]. We explore the 
normalized current skewness F3 = K3/K1 versus both the ballistic-diffusive cross-over and the 
thermal (fc^T 3> eV) - shot noise (eV 3> fcgT) cross-over. At zero temperature, it was shown[42] 
that N as small as 8 gives a correct picture of the skewness, even in the diffusive limit. We 
performed a numerical simulation of such a small N by solving numerically the largest eigenvalue 
of the counting matrix. Fig. 9 shows F^{L / L e i as t, eV/ksT). For comparison, the continuous line is 
calculated from the analytical results in the diffusive limit and good agreement is found with the 
numerical estimation. This is an interesting result, which shows a tendency of the simple exclusion 
model to reproduce the same FCS as for a larger class of models. This raises the question of trying 
to identify more precisely the key features which are required for a stochastic model to yield this 
FCS. The pearly line corresponds to a single scatterer with the same conductance as the system 
in the low temperature limit. As L/L e i ast — > {L may be arbitrary large), the system is expected 
to converge to such a point scatterer limit and a good agreement is found with the simulation. 
Good agreement is also found with Monte-Carlo simulations performed along the ballistic-diffusive 
cross-over at zero temperature [42]. 



4 Connection to microscopic models 

As we have seen in section 3 above, there are many instances where a simple approach in terms of 
purely classical exclusion models yields an impressive agreement with fully quantum-mechanical 
computations, specially for diffusive systems. This raises naturally the question of trying to 
derive these classical models directly from a fully quantum mechanical description of a coherent 
conductor. In this section, we shall establish such a connection for the counter-flows exclusion 
model, but only in the rather special situation of a single barrier. 



The basic principle of this model is to discretize the single-particle states of conduction elec- 
trons. For each value of the kinetic energy, we consider two types of localized wave-packets, moving 
in either direction (from right to left, or left to right). Each possible wave-packet is supposed to be 
spatially located in one of N cells of width Ax into which our mesoscopic conductor is subdivided. 
Let us first assume a perfect conductor, with no impurities. It is then natural to approximate the 
true quantum mechanical evolution operator describing a single electron during a time interval 
At = Ax/vp (vp is the Fermi velocity of electrons, taken to be also the group velocity of these 
wave-packets) by the following matrix U acting as: 

U\i,+) = \i + l,+) (42) 
U\i,-) = (43) 

Here, \i, ±) denotes a wave-packet with velocity ±vp located in the cell labelled by the integer i. 
From now on, we shall assume i ranges over all integers, ignoring what happens at the boundaries 
of the sample. In fact, rather than taking the reservoirs explicitely into account here, we have in 
mind an evolution of the system in a finite time interval, from an initial state with a finite number 
of particles, located on the left and on the right of the finite length region where impurity scattering 
may occur. After all, reservoirs always contain a finite number of particles in real experiments. It is 
straightforward to check that U thus defined is unitary. States with p electrons located in cells i\ < 
i 2 < •■• < i p and moving with velocities v\, V2, v p (vi = ±) are obtained as usual by forming fully 
antisymmetrized tensor products of the single particle states \ii,Vi), |i 2 j ^2), \i P , v p ), denoted by 
|*i , vi; i2,V2', --^ip, v p ). The unitary operator U allows to construct the corresponding evolution 
operator in the space of p electron states, and this latter operator is defined by: 

U {p) \i 1 ,v 1 ;i 2 ,v 2 ; ...;ip,v p ) = \h + v 1 ,v 1 ;i 2 + v 2 ,v 2 ; ...;i p + v p ,v p ) (44) 

Now, if we prepare such a system at time iina density matrix p(t) which is diagonal in the above- 
given basis of p electron states, we have p(t + At) = p(t)U^ + , and equation( 44) shows that p 
remains diagonal. This defines a very simple process, where the probability V(ii, v\; 12, v 2 ; ...;ip,v p ;t) 
evolves according to: 

V(ii,v 1 ;i 2 , v 2 ; ...; i p , v p ; t + 1) = V(h -v 1 ,v 1 ;i 2 - v 2 , v 2 ; i p - v p , v p ; t) (45) 
Note that time is now written in units of At, which simplifies the notation. 

The next step is to include a single barrier, located between cells i and i + The definition of 
the single particle evolution operator U is now modified according to: 

U\j,+) = \j + 1,+) (j^i-1) (46) 

U\j,~) = (j^i) (47) 

U\i-l,+) = cos(0i)|i,+)+sin(0j)|i-l,-) (48) 

U\i,-) = cos(6i)\i-l,-)-sm(6i)\i,+) (49) 

There is now a finite probability sin 2 (#i) for the particle coming from the left (on cell i — 1) or 
from the right (on cell i) of the barrier to be backscattered. Let us denote by the transmission 
probability, equal to cos 2 (#i). Note that the U matrix just written is not the most general unitary 
matrix, and both transmission and reflection amplitudes are in general complex numbers. But 
again, we made this choice to keep a lighter notation. As before, this model for the single particle 
quantum evolution generates also the evolution of p particle states. Since the only single particle 
states which are affected by the presence of this barrier are \i — 1, +) and \i, — ), it is just sufficient 
to consider initial states where these are both occupied. In fact, because of the Pauli principle, 
we have: 

UW\i-l,+;i,-) = \i,+;i-l,-) (50) 

We shall not write explicitely the corresponding unitary matrix here. But the important point 
to stress is that it induces an evolution of the density matrix p which still preserves the stability of 
the subspace of diagonal matrices. The basic idea behind this property is a description of the time 
evolution of the system with p particles in terms of trajectories, in the spirit of the path integral 
approach to quantum mechanics. In the presence of a unique barrier, there is at most a single 
path which connects an initial state t>i; 12, V2; ip, v p ) to a final state v[; i 2 , v 2 ; i' p , v' p ). 
This holds because under successive applications of U, the distance of a particle to the barrier in 
the initial state \i,v) which evolves into a known final state \i',v'), is uniquely determined by this 
final state and the total duration of this evolution. So if the distances of the final locations of the 
particles ii,i' 2 , ...,i' p to the barrier are all different, there is only one permutation of the particles 
which takes the i's into the i"s. When two i n s correspond to the same distance to the barrier, 
the Pauli principle imposes these two sites to be located on opposite sides of the barrier, and the 
same statement is true for the i's of the corresponding initial states. If we had distinguishable 
particles, this would yield two possible paths connecting this pair of initial i's to this pair of final 



i"s. But the antisymmetrization operation forces us to see these two paths as a single one, as in 
equation (50) above. 

As for the pure system, we may then interpret the evolution of the diagonal part of the density 
matrix in terms of a simple stochastic process for the probability V(i\, v\; 12, V2', ■ i p , v p ; t), which 
is nothing but the counter-flows exclusion model with a single barrier, defined in section 2 above. 

Unfortunately, it is difficult to generalize this construction to more than a single barrier. For 
three barriers and a single particle already, there are in general several differents paths (generated 
by a given number of successive applications of the matrix U) which connect an initial state \i,v) 
to a final state \i',v'). Therefore, the corresponding evolution of the density matrix no longer 
preserves the subspace of diagonal matrices. 

We took here the view that we may still obtain some valuable information on the FCS by 
ignoring the non-diagonal elements of the density matrix. This approximation is rather severe for 
a long coherent single channel conductor with a large number of barriers, since it does not allow 
to capture the phenomenon of Anderson localization. But the good agreement obtained between 
purely classical models and fully quantum ones in the case of a two barriers and of a diffusive 
medium, certainly calls for a deeper understanding of the possible connections between these two 
classes of models. 



5 Conclusion 

In this paper, we have checked explicitely that the FCS of purely classical stochastic models 
with a local exclusion constraint coincides with the FCS of a large class of quantum mechanical 
microscopic models where the current flows between two reservoirs. In particular, no restriction 
on the system size seems to be necessary, since this coincidence appears for a small number of 
tunnel barriers, as well as in the limit of a diffusive medium composed of an arbitrary large number 
of such barriers. The key common ingredient to both types of models is the constraint removing 
doubly occupied sites, imposed by the Pauli principle. 

It would be very interesting to extend the classical approach in terms of exclusion models 
to multi-terminal geometries. Such generalizations have been considered in a fully quantum- 
mechanical treatment of diffusive systems [6], or in the semi-classical Boltzmann-Langevin ap- 
proach [47]. These works have so far focused on the correlation matrix of the integrated charges 
flowing through each contact during a finite time interval. A natural question is whether simple 
classical exclusion models are able to generate these results for the second order cumulants. The 
computation of higher-order cumulants in multi-terminal geometries is another interesting open 
question. 
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